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We study the sliding friction for viscoelastic solids, e.g., rubber, on hard flat substrate surfaces. 
We consider first the fluctuating shear stress inside a viscoelastic solid which result from the ther- 
mal motion of the atoms or molecules in the solid. At the nanoscale the thermal fluctuations are 
very strong and give rise to stress fluctuations in the MPa-range, which is similar to the depinning 
stresses which typically occur at solid-rubber interfaces, indicating the crucial importance of ther- 
mal fluctuations for rubber friction on smooth surfaces. We develop a detailed model which takes 
into account the influence of thermal fluctuations on the depinning of small contact patches (stress 
domains) at the rubber-substrate interface. The theory predict that the velocity dependence of the 
macroscopic shear stress has a bell-shaped form, and that the low- velocity side exhibit the same 
temperature dependence as the bulk viscoelastic modulus, in qualitative agreement with experimen- 
tal data. Finally, we discuss the influence of small- amplitude substrate roughness on rubber sliding 
friction. 



1. Introduction 

The friction between rubber and smooth substrate sur- 
faces is a topic of extreme practical importance, e.g., for 
wiper blades (in particular on hydrophobic glass), rub- 
ber 0-ring seals, and in the contact region between the 
tire-rubber and the steel rim on a wheel Q. 

When a rubber block is sliding on a very rough sub- 
strate, such as a tire on a road surface, the friction is 
almost entirely due to the energy dissipation in the bulk 
of the rubber as a result of the fluctuating (in time and 
space) viscoelastic deformations of the rubber by the sub- 
strate asperities 2, 3, 4J. This mechanism becomes unim- 
portant when the substrate is very smooth. In the limit- 
ing case of a perfectly smooth substrate, the friction is in- 
stead due to local stick-slip events at the sliding interface. 
Schallamach0 has proposed a molecular mechanism for 
the local stick slip, see Fig. ^ where rubber polymer 
chains at the interface attach to the moving countersur- 
face, stretches, detaches, relaxes, and reattaches to the 
surface to repeat the cycle (similar models have been 
studied in Ref. IE 0)- During each cycle, the elastic 
energy stored in the polymer chain is dissipated as heat 
during the (rapid) detachment and relaxation phase, and 
this is assumed to be the origin of the (macroscopic) fric- 
tion. However, in our opinion this picture cannot be 
fully correct. First, for physisorption systems the en- 
ergy barriers for (vertical) detachment are usually much 
higher than the energy barriers for lateral sliding so 
one would not expect any detachment to occur. Secondly, 
with respect to the stresses rubber materials are usually 
exposed to, rubber is nearly incompressible and it is not 
easy to imagine how single molecules strongly confined 
at the interface are able to switch between an elongated 
(stretched) state and a relaxed (curled up) state as indi- 
cated in the figure. Furthermore, the sliding friction tend 
to exhibit the same temperature dependence as the bulk 
rubber viscoelastic modulus E(uj) as described by the 
Williams-Landel-Ferry (WLF) 9] shift factor ar- This 
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FIG. 1: The classical description of a polymer chain in 
contact with a lateral moving count ersur face. The chain 
stretches, detaches, relaxes, and reattaches to the surface to 
repeat the cycle. Adapted from Ref. |20|. 



indicate the crucial role of the bulk rubber in the friction 
process. 

We believe that the local stick-slip processes, which 
must occur at the sliding interface, involve relative larger 
rubber volume elements, always in adhesive contact with 
the substrate. That is, during sliding small patches 
(with a diameter of order D ^ 10 — 100 nm) or stress 
domains 10] of rubber at the sliding interface perform 
stick-slip motion: during stick the shear stress at the in- 
terface increases continuously with time until the local 
shear stress reaches a critical depinning stress dc, after 
which a rapid local slip occur, but with the rubber patch 
in continuous adhesive contact with the substrate. Dur- 
ing the local slip the elastic deformation energy stored 
in the rubber during the loading phase will be dissipated 
party inside the rubber (in a volume of order D^) and 
partly at the interface. The deformation field in the 
vicinity of a stress domain of area ~ will extend a 
distance ~ D into the rubber block; we denote this basic 
unit (volume ~ D^) as a stress block [see volume elements 
surrounded by dashed lines in Fig. I2ta)]. 

Figure [21 illustrate three other mechanisms of rubber 
friction, which all depend on the rubber-substrate ad- 
hesional interaction. Fig. Eljb) illustrate a case where a 
rubber block is sliding on a smooth wavy substrate. Here 
we assume only roughness (waviness) on a single length 
scale, i.e., the substrate bumps have no roughness on 
length scales smaller than the lateral size of the bumps. 
In this case, if the adhesional interaction (and the exter- 
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FIG. 2: For different processes where the adhesive rubber- 
substrate interaction contribute to rubber friction. 



nal applied normal load) is unable to bring the rubber 
into perfect contact with the substrate, (as in the fig- 
ure), during sliding opening and closing cracks will occur 
at the exit and the front of each asperity contact region. 
In general, negligible bulk viscoelastic energy dissipation 
occur at the closing crack, while a huge energy dissipa- 
tion may occur at the opening crack[ll[ [ij. It has been 
shown that in some cases the dominant contribution to 
the friction force arises from the energy "dissipation" at 
the opening 

cracks [illE3- 
If the substrate has small-amplitude roughness, the 
adhesive rubber-substrate interaction may (even in the 
absence of an external load) lead to complete contact be- 
tween the rubber and the substrate at the sliding inter- 
face [see Fig. Etc)]. In this case, during sliding viscoelas- 



FIG. 3: The thermal motion of the atoms in a solid gives rise 
to shape fluctuations for a small volume element in the solid 
(schematic). 



tic deformations will occur in the bulk of the rubber in 
the vicinity of the substrate, which will contribute to the 
observed friction. This contribution can be calculated 
using the theory developed in Ref. 0, E| • 

Finally, for very soft rubber it has often been observed 
that detachment waves propagate throughout the entire 
contact area, from its advancing to the trailing edge, see 
Fig. Eljd). Schallamach gl4j first discovered these waves at 
"high" sliding velocities and for (elastically) soft rubber. 
Roberts and Thomas [HI [l5| have shown that when such 
instabilities occur, the frictional stress is mainly due to 
the energy dissipation at the opening crack, i.e., similar 
to the situation in Fig. Eljb). 

In this paper we study the rubber friction process 
shown in Fig. Efa). This is probably the most important 
rubber friction mechanism in most applications involving 
very smooth surfaces. In Sec. 2 we consider the fiuctu- 
ating shear stress inside a viscoelastic solid, which result 
from the thermal motion of the atoms or molecules in 
the solid. At the nanoscale the thermal fiuctuations are 
very strong giving stress fiuctuations in the MPa-range, 
which is similar to the depinning stresses which typically 
occur at solid-rubber interfaces, illustrating the crucial 
importance of thermal fiuctuations for rubber friction on 
smooth surfaces. In Sec. 3 we develop a detailed model 
which takes into account the infiuence of thermal fiuc- 
tuations on the depinning of small contact areas (stress 
domains) at the rubber-substrate interface. The theory 
predict that the velocity dependence of the macroscopic 
shear stress has a bell-shaped form (see Sec. 4), and that 
the low-velocity side exhibit the same temperature de- 
pendence as the bulk viscoelastic modulus, in qualitative 
agreement with experimental data (Sec. 5). In Sec. 6 
we discuss the role of the surface roughness, which exist 
even on very smooth surfaces. We show that the coun- 
tersurfaces used in most rubber sealing applications have 
such large roughness that the friction observed during 
sliding at low velocity in many cases may be mainly due 
to the substrate asperity-induced viscoelastic deforma- 
tion of the rubber surface. Sec. 7 contains the summary 
and conclusion. 



2. Brownian motion in viscoelastic solids 

A small particle in a liquid perform random motion 
caused by the impact of the surrounding liquid molecules. 
Similarly, a small volume element in a viscoelastic solid 
perform shape fluctuations as a result of the thermal mo- 
tion of the atoms or molecules in the solid, see Fig. [HI 
Here will will estimate the magnitude of the fluctuating 
shear stress which act on any (internal) surface in a vis- 
coelastic solid. 

Assume that u(x, t) is the displacement vector in an 
infinite viscoelastic solid. The equation of motion for u 
is 

=AV'u+(A + A)VV-u + f (1) 

where fi and A are time integral operators, e.g., 

fiG{t) = j dt' ii{t-t')G{t'), 

and where f (x, t) is a randomly fiuctuating force density 
(see below). If we define the Fourier transform 
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Let us study the fiuctuating stress in the solid. The aver- 
age stress {aij ) = but the average of the square of any 
components of a^j will be finite. Here we consider the 
magnitude of the fiuctuating shear stress within some 
plane which we can take to be the xy-plane. Thus we 
consider 
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Using (2) and (3) and assuming that v is frequency in- 
dependent gives 
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In deriving (4) we have used the relation ji = E/2{l^v)^ 
where E{(jo) is the Young's (viscoelastic) modulus. The 
shear stress when we only include wavevectors up to = 
2t:/D (where I) is a cut-off length, of order the distance 
between the cross-links, or larger) is given by 



= -k^TCD- 
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where Eq = E{0) and = ^(oo) are the low- and 
high-frequency rubber modulus (both real), we get 

^\\ = H) = k^TCD-\E^ - Eo) 

Assuming u = 0.5 gives 27rC/3 ~ 1.1. In a typical case 
Eoo ^ 10^ Pa >> £^0 and if i:> = 10 nm we get at 
room temperature ay ^ 1 MPa, which (typically) is of 
the same order of magnitude as the depinning stress at a 
rubber-substrate interface. 



It is also interesting to estimate the fluctuation in the 
displacement uy and the fluctuation in the strain. Using 
the same approach as above one obtain 
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we get 
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In a similar way one can calculate the average of the 
square of the strain 

{i) = {el, + e%) = k^TC"D-' - ^ 
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With Eq = 1 MPa and D = 10 nm one get at room 
temperature u\\ ~ 1 nm and ey ^ 0.3. 

The analysis above indicate that, on the length scale of 
~ 10 nm, very large strain and stress fluctuations occur 
in normal rubber. The stress fluctuations are of similar 
magnitude as the (typical) rubber-substrate depinning 
stress, suggesting that thermally excited transitions over 
the (lateral) energy barriers play a crucial role in rub- 
ber friction on smooth substrates. However, the stress 
fluctuations (~ 1 MPa) are negligible compared to the 
stress necessary for detaching a rubber patch in the nor- 
mal direction, since the latter is determined by the ad- 
hesional stress 0,^3 which typically is of order 1 GPa. 
Thus, rubber sliding on smooth surfaces never involves 
(thermally activated) detachment of rubber from the sub- 
strate surface, but only (lateral) sliding of rubber patches 
(the lateral depinning stress is typically 100 — 1000 times 
smaller than the adhesional stress). [Detached areas may 
form and propagate at the interface (as for Schallamach 
waves, see flg. EJd)), but in these cases they are gener- 
ated by the external applied stress, and depend on the 
shape of the bodies; the detached regions usually form at 
the edge of the rubber-substrate contact region, where 
elastic instabilities (e.g., "buckling") of the rubber may 
occur 18j.] In the next section we will develop a model 
of rubber sliding friction based on the picture presented 
above. 



3. Theory 

We consider a rubber block with a smooth flat surface 
sliding on a perfectly smooth and flat substrate. We as- 
sume that the adhesive rubber-substrate interaction will 
result in perfect contact between the two solids (i.e., we 
assume that no Schallamach waves occur at the inter- 
face). During sliding at low velocities, small (nanometer 
sized) regions (stress blocks) at the sliding interface will 
perform stick-slip motion. We model the real physical 
system in Fig. ^a) with the block-spring model shown 
in Fig. ^b). However, the springs in the model are 
not elastic springs but viscoelastic springs determined by 
the viscoelastic modulus of the rubber (see below). Fur- 
thermore, the blocks in the model experience not only 
the stress from the substrate and the forces from the 
springs, but also randomly fluctuating (in time) forces 
derived from the thermal motion of the molecules in the 
solid. The strength of the fluctuating forces is determined 
by the temperature and by the viscoelastic properties of 
the springs via the fluctuation-dissipation theorem. The 
combination of the (thermal) fluctuating forces, and the 
forces derived from the external pulling of the upper sur- 
face of the rubber block, and the stress acting on the 
bottom surface of the block from the substrate, deter- 
mine the motion of the stress blocks. 

3.1. Basic equations 

We assume that the motion at low sliding velocity oc- 
cur by a thermally activated process, where small sur- 
face areas or "patches" of rubber at the interface per- 
form stick-slip motion. We will refer to these "patches" 
as stress domains. The stress domains are pinned by the 
substrate potential, and we assume that some character- 
istic shear stress dc must be reached before local slip can 
occur. Thus, the pinning force Fc = (TcD^^ where D is 
the characteristic linear size of a stress domain. 

The equation of motion for the coordinate qi of the 
stress domain i is assumed to be 

mqi = kz{x- Qi) + kx{qi+i - Qi) 



+^x(<7i-i - Qi) ^ fi^Fi 
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where x = vt and where kz is a time integral operator 
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and similar for kx- We consider N stress domains, 
i = 1,...,7V, and assume periodic boundary conditions 
so that qo = Qn- In (5), fi{t) is a stochastically fluctu- 
ating force which result from the thermal motion of the 
rubber molecules. The force Fi act on the stress domain 
i from the substrate and is deflned as follows: When the 
stress domain i slips, then Fi = —mrjqi. In the pinned 
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FIG. 4: (a) Rubber block in adhesive contact with a flat sub- 
strate. During sliding at low velocities, small rubber volume 
elements (stress blocks) perform stick-slip motion, (b) The 
model used in the mathematical description of the system in 
(a). The viscoelastic springs k are determined by the rubber 
viscoelastic modulus and by the lateral size D of the stress 
blocks via k((jj) ~ E{ijj)D. 



state, Fi is just large enough to balance the total force 
exerted by the stress block on the substrate: 



-kz{x - Qi) - kx{qi^i - qi) - kx{q, 



'i-i 



Slip will start when \Fi \ reaches a critical value Fc, either 
as a result of the external applied force or as a result of 
a large enough thermal fluctuation, or, in general, as a 
combination of both these effects. 
We define the Fourier transform 



x{uj) = — [ dt x{t)e'^' 
2ti J 

x{t) = [ duj x{uj)e-''^' 



(8) 
(9) 



The fluctuating force fi(t) result from the thermal motion 
of the rubber molecules and must satisfy the fluctuation- 
dissipation theorem. That is, if we write 

Kijqj = kzqi - kx{qi+i - qi) - kx{qi-i - qi) (10) 

then from the theory of Brownian motion (see Appendix 
A) 

{fi{i0)fj{L0')) = -^luiKij{L0)5{L0 + LO') (11) 
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and (for t > 0): 

Kij{t) = ^ J duj Kij{uj)e~ 
If we write the elastic modulus asl28l 
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Eq. (12) gives 
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Substituting (13) in (11) gives 
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We can write 

-Kij{u)qj{uo) + kz{uo)x{uo) 

= -K^qjico) + k^xico) -^hnUmico) (15) 
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Here we have added a stochastically fluctuating force on 
the right hand side of (17) which we can choose so as to 
reproduce the fluctuating force fi{t) in Eq. (5). That 
is, if we choose Qniif) appropriately, we can remove the 
force fi{t) in Eq. (5). To this end we must choose 

fM = -Y.T^ (18) 
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Using (18) and (27) gives 
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which agree with (14). 

Let us summarize the basic equations: 

mqi = k^{x- Qi) + k^iqi^i + Qi-i - 2qi) 



where 
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are complex Gaussian random variables: where 
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Using (19) and (22) gives 
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Using (23)- (25) this gives 

{gnl{(^)gn'l'{(^')) 

= A^~"^ M^^cos[A:(x/ — xif)]Snn'^{^ + where 



3.2. Numerical implementation 

If D is the lateral size of a stress block, then the mass 
of a stress block m = pD^. We introduce the spring 
constants /c^ = a^k* and k^ = axk* where k* = D£^oo, 
and where and ce;^ are of order unity. We measure 
time in units of r = (m//c*)-^/^ and distance in unit of 
/ = Fc/k* where Fc = cFcD'^ is the stress block pinning 
force. We also measure u in units of Fc and x in units of 
/. In these units we get 
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In order to numerically integrate the equations (33) and 
(34) we discretize time with the step length 5t. The (fluc- 
tuating) force 6gni to be used for each time-step in (34) 
can be written as 



But if 
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we get 
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Thus we can write = (27r^t)^/^G where G is a Gaus- 
sian random number with (G^) = 1. Thus, we take 



where G^^\ and G^^^ are Gaussian random numbers. We 
can also write 

^9ni = N-y^Re Mfcne"=^^ (G^ + iCfj ) (35) 
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In the calculations presented below we have assumed 
that a sliding block return to the pinned state when the 
shear stress < ad = XcTc^ where cic = Fc/D^ is the 
depinning stress and A < 1. We assume that when the 
shear stress has decreased below (Jd the block return to 
the pinned state with the probability rate and we use 
random number to determine when the transition actu- 
ally takes place. Thus, if r is a random number uniformly 
distributed in the interval [0, 1], then if \a\ < ad we as- 
sume that the stress block return to the pinned state dur- 
ing the time interval St (the time integration step length) 
if wSt > r. In the simulations below we use w < 10^^ s~^ 
and A = 0.1. We use 5 ~ 10~^^ s in most of our simu- 
lations so that the condition wSt << 1 is satisfied, and 
the results presented below does not depend on the time 
step 5t. The calculation of Sgnj (Eq. (35)) is conveniently 
performed using the Fast Fourier Transform method. 

4. Numerical results 

We now present numerical results for styrene- 
butadiene rubber both with and without filler. 

Fig. Efa) shows the frictional shear stress df and (b) 
the fraction of slipping surface area N/Nq as a function 
of the logarithm (with 10 as basis) of the sliding velocity 
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FIG. 5: (a) The frictional shear stress fjf and (b) the fraction 
of slipping surface area N/Nq as a function of the logarithm 
(with 10 as basis) of the sliding velocity of the rubber block. 
For Styrene butadiene copolymer with 60% carbon black and 
for two different temperatures, T = 10 °C and 20 °C. In the 
calculation we have used the (zero-temperature) depinning 
stress cTc = 1 MPa and the stress block size D — 30 nm. The 
number of stress blocks was 128. The viscous friction coeffi- 
cient during steady sliding 77 = 0.03 (natural units) and the 
critical stress below which steady sliding becomes metastable 
(Jci = O.lcTc =0.1 MPa. The probability rate per unit time to 
return to the pinned state when a < ad is w = 2 x 10"^° s~"^. 



of the rubber block. The results are for styrene butadiene 
copolymer with 60% carbon black, and for two different 
temperatures, T = 10 °C and 20 °C. Note that the low 
velocity part {v < vi, where vi ^ 10~^ m/s is the velocity 
at which the friction is maximal) of the friction curve is 
shifted by ~ 1 decade towards lower velocities when the 
temperature is reduced from 20 to 10 °C. This is iden- 
tical to the change in the (bulk) viscoelastic shift factor 
ttT, which change by a factor of ~ 10 during the same 
temperature change. Fig. ]^h) shows that the number 
of moving stress blocks is basically temperature indepen- 
dent for sliding velocities v > vi ^ 10~^ m/s, i.e., for 
V > vi the fluctuating force arising from the finite tem- 
perature has a negligible influence on the friction force. 
For V < vi more stress blocks are depinned at the higher 
temperature because the fluctuating force is larger at the 
higher temperature. Note also that the maximal shear 
stress cTf ^ 0.3 MPa is about 30% of the depinning stress 
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FIG. 6: The frictional shear stress df as a function of the 
logarithm of the shding velocity of the rubber block. For the 
temperature T = 20 °C and for the stress block sizes D = 30 
and 40 nm. All other parameters as in Fig. 5. 



(7c = I MPa. 

Fig. inija) shows the frictional shear stress df and (b) 
the fraction of sliding blocks, as a function of the loga- 
rithm of the sliding velocity. In the calculation we have 
used T = 20 °C and the stress block sizes I) = 30 and 
40 nm. When the stress block size increases the pin- 
ning force cFcD'^ increases, and it is necessary to go to 
lower sliding velocities in order for temperature effects to 
manifest them-self (as a decrease in the frictional shear 
stress). 

The linear size of the stress blocks, is most likely de- 
termined by the elastic modulus and the depinning shear 
stress (Jc as follows. The stress block is the smallest unit 
which is able to slide as a coherent unit and can be de- 
termined as follows. The depinning force is Fc = cFcD^. 
If the stress a act at the bottom surface of a stress block 
it will move a distance x determined by kx = crD^ where 
k ~ ED^ where E is the elastic modulus. The stress 
block experience a quasi-periodic pinning potential from 
the substrate, characterized by a lattice constant a of 
order a few Angstroms. Thus, the stress block will in 
general be able to occupy a "good" binding position in 
the corrugated substrate potential only if ka is less than 
the pinning potential The condition ka = Fc gives 
the size D of the pinned domains. Using Fc = cFcD^ and 
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FIG. 7: The frictional shear stress (jf as a function of the 
logarithm of the sliding velocity of the rubber block. For the 
temperature T = 20 °C and for the viscous friction coefficients 
rj = 0.03 and 0.01 (natural units). All other parameters as in 
Fig. 5. 



k ~ ED this gives 



D ^ Ea/cFc. 



In a typical case E ^ 1{) MPa, <Jc ~ 1 MPa and a a 
few Angstroms, giving D = 30 nm. Since D increases 
when the elastic modulus increases, the present theory 
indicate that the friction should increase with increasing 
elastic modulus. This is exactly what has been observed 
in friction studies for smooth surfaces (see below). Since 
the elastic modulus E depends on frequency, during slid- 
ing the size D of the stress domains may depend on the 
sliding velocity, but we have not taken this effect into 
account in this paper. 

Fig. [3 shows the frictional shear stress df as a function 
of the logarithm of the sliding velocity. Results are shown 
for two different viscous friction coefficients t] = 0.03 and 
0.01 (natural units). As expected, rj is only important at 
relative high sliding velocities. 

Fig. IHl shows how df depend on the sliding velocity 
for two different values (0.1 and 0.2 MPa) of the critical 
stress CTci below which the sliding patch can return to the 
pinned state. 

Fig. ini shows the frictional shear stress af as a function 
of the logarithm of the sliding velocity of the rubber block 
for two different values (2 x 10^^ and 10^^ s~^) of the 
probability rate per unit time to return to the pinned 
state when a < <Jci. Fig. ITOl shows similar results for two 
different values (1 and 0.7 MPa) of the depinning stress 

CTc 

Fig^J shows the frictional shear stress df as a function 
of the logarithm of the sliding velocity of the rubber block 
for an unfilled styrene butadiene (SB) copolymer for two 
different values (30 and 40 nm) of the size D of a pinned 
region. 
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FIG. 8: The frictional shear stress at as a function of the loga- 
rithm of the shding velocity of the rubber block. For the tem- 
perature T = 20 °C and for the critical stresses (below which 
steady sliding becomes met ast able) ad = 0.1 ctc = 0.1 MPa 
and 0.2 MPa. All other parameters as in Fig. 5. 



5. Discussion 

Vorvolakos and Chaudhury 20j have performed a very 
detailed experimental study of sliding friction for silicon 
rubber on hard flat (passivated) substrates (see also Ref. 
hi , 22, 23] for other studies of elastomer friction). They 
used silicon rubbers with many different (low-frequency) 
elastic modulus E. 

In Fig. El we show the velocity dependence of the 
shear stress as measured at different temperatures T = 
298 (open circles), 318 (gray circles) and 348 K (black 
circles) for a silicon rubber with the low-frequency elastic 
modulus E ^ 5 MPa. The experimental data on the 
low-velocity side can be shifted to a si ngle curve when 
plotted as a function of var^ see Ref. I2a. This IS m 
accordance with our model calculations (see Sec. 4) and 
shows the direct involvement of the rubber bulk in the 
sliding dynamics. 

Vorvolakos and Chaudhury 20j have shown that the 
frictional shear stress for all the studied rubbers increases 
with increasing elastic modulus. This is in qualitative 
agreement with our theory since, as explained in Sec. 4, 
as E increases we expect the linear size D of the stress 
domains to increase, which will increase the sliding fric- 
tion (see Fig. 151 and [TT|) . 

In the experiments by Groschji^ it was observed that 
the friction on smooth surfaces has a bell-like shape with 
a maximum at some characteristic velocity v = vi. We 
observe the same general behavior, but with some impor- 
tant differences. Thus, the experimental data of Grosch 
was observed to obey the WLF transform. That is, the 
full friction curve could be constructed by performing 
measurements in a very limited velocity range and at dif- 
ferent temperatures, and then use the WLF transform to 
shift the measured data to a single temperature. We also 
find that our calculated friction obey the WLF transform 
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FIG. 9: The frictional shear stress df as a function of the 
logarithm of the sliding velocity of the rubber block. For 
the temperature T = 20 °C and for the probability rates 
per unit time to return to the pinned state (when a < ad) 
w — 2 X 10"*^° s""*^ and w — 10"*^^ s""*^. All other parameters as 
in Fig. 5. 




-6 -4 -2 

log V (m/s) 



FIG. 10: The frictional shear stress df as a function of the 
logarithm of the sliding velocity of the rubber block. For 
the temperature T = 20 °C and for the (zero-temperature) 
depinning stresses ctc = 1 MPa and and ac = 0.7 MPa. All 
other parameters as in Fig. 5. 
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FIG. 11: The frictional shear stress fjf as a function of the 
logarithm of the shding velocity of the rubber block. For un- 
filled styrene butadiene (SB) copolymer for T = 20 °C. In 
the calculation we have used the (zero-temperature) depin- 
ning stress ctc = 1 MPa and the stress block size D = 30 nm 
and 40 nm. The number of stress blocks was 128. The vis- 
cous friction coefficient during steady sliding 77 = 0.03 (natu- 
ral units), and the critical stress below which steady sliding 
becomes metastable ad = O.Ictc = 0.1 MPa. The probability 
rate per unit time to return to the pinned state when a < ad 
isw = 10^2 s-^ 
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FIG. 13: The frictional shear stress df as a function of the 
logarithm of the sliding velocity of the rubber block. For 
Styrene butadiene copolymer with 60% carbon black and for 
T = 40 °C. The friction is entirely due to the surface rough- 
ness of the substrate which is assumed to be self affine frac- 
tal with the fractal dimension Df — 2.3 and the root-mean- 
square roughness 0.5 nm. The long- wavelength and short 
wavelength roll-off wavevectors qo = 27r/Ao and qi = 27r/Ai 
where Ao = 100 /xm and Ai = 2 nm (upper curve) and 10 nm 
(lower curve). 
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FIG. 12: Shear stress as a function of velocity and temper- 
ature for a silicon elastomer (low frequency elastic modulus 
E ^ 5 MPa) sliding on a Si wafer covered by an inert self- 
assembled monolayers film. Open circles, gray circles and 
black circles represent data at 298, 318 and 348 K, respec- 
tively. Adapted from Ref. 

for V < vi, but not for v > vi. Thus, the decrease in the 
friction which we observe for v > vi is nearly temper- 
ature independent. The difference between our predic- 
tion and the Grosch results can be understood as follows: 
Most likely, the decrease of the friction for v > vi in the 
Grosch experiment is related to a decrease in the area of 
real contact with increasing sliding velocity. The Grosch 



experiments where performed on smooth but wavy glass, 
and the area of real contact depends on the effective elas- 
tic modulus of the rubber. Thus, as the sliding velocity 
increases or, equivalent ly, the temperature decreases, the 
rubber becomes more stiff and the area of real contact 
decreases. Hence, if the shear stress remains constant at 
high sliding velocity [as we indeed observe in our calcu- 
lations if the calculations are performed at low velocities 
and different temperatures, and shifted to higher veloc- 
ity using the WLF equation], then the friction force will 
decrease with increasing sliding velocity because of the 
decrease in the area of real contact. 

6. The role of surface roughness: application to 
rubber sealing 

The experiments by Vorvolakos and Chaudhurv|2oj 
was performed on (passivated) silicon wafers with the 
root-mean-square roughness of at most 0.5 nm. For such 
smooth surfaces the adhesive rubber-substrate interac- 
tion will, even in the absence of a squeezing pressure, 
give rise to complete contact between the rubber and the 
substrate within the nominal contact area. For this case 
we have applied the rubber friction theory developed in 
Ref. to obtain the contribution to the frictional shear 
stress from the roughness induced bulk viscoelastic de- 
formations in the rubber (Fig. [2h). We have assumed 
that the substrate is self affine fractal with the fractal 
dimension = 2.3 and the long-distance roll off wave- 
length Ao = 100 fim. We have included roughness compo- 
nents down to the short-distance (cut-off) wavelength Ai. 




-6 -4 -2 2 4 
log V (m/s) 



FIG. 14: The kinetic friction coefficient /ik as a function of 
the logarithm (with 10 as basis) of the shding velocity of the 
rubber block. For Styrene butadiene copolymer with 60% 
carbon black and for T = 40 °C. The friction is entirely due 
to the surface roughness of the substrate which is assumed to 
be self affine fractal with the fractal dimension Df = 2.3 and 
the root-mean-square roughness 1 /xm. The long-wavelength 
and short wavelength roll-off wavevectors = Stt/Aq and 
qi — 27r/Ai where Ao — 100 /xm and Ai = 2 nm (top curve), 
10 nm (middle curve) and 50 nm (bottom curve). The solid 
curves are the full calculation while the dashed curves are 
without the flash temperature effect. 

In Fig. ^|we show the resulting frictional shear stress 
(assuming perfect contact at the sliding interface) using 
Ai = 10 nm and 2 nm. The smallest possible Ai is de- 
termined by an atomic distance but a more likely cut-off 
length in the present case is the mean distance between 
rubber cross links, which is of order a few nanometer. 
The calculated maximal shear stress is of order ~ 100 Pa, 
which is a factor ^10^ smaller than the shear stress ob- 
served in Ref. ^2^. Hence we conclude that the surface 
roughness in the measurements of Vorvolakos and Chaud- 
hury has a negligible influence on the observed frictional 
shear stress. However, the surfaces used in Ref. 20j are 
exceptionally smooth, and it is of great interest to study 
the influence of the surface roughness on rubber friction 
for surfaces of more common use in rubber applications. 

One very important application of rubber is for seals. 
In a typical such application the rubber is sliding on a (lu- 
bricated) steel surface with the root-mean-square rough- 
ness of order 1 /im. At low sliding velocity (e.g., at the 
start of sliding) the lubricant fluid is squeezed away from 
the rubber-substrate asperity contact regions, giving fric- 
tion coefficients similar to those of dry (unlubricated) 
surfaces, typically of order /ik ~ 0.5. 

Recently Mofidi et alji^l have performed rubber fric- 
tion studies for nitrile rubber (the most popular seal ma- 
terial) in contact with a steel surface. The experiments 
was performed by squeezing a steel cylinder (diameter 
1.5 cm and length 2.2 cm) with the force = 100 N 




FIG. 15: Rubber in squeezed contact with a hard rough 
substrate. During sliding the bulk viscoelastic deformation 
(on the length scale A) at an asperity contact region will not 
be influenced by the thin lubrication film (thickness h) at the 
asperity if /i << A. The lubrication film will, however, remove 
the contribution to the friction from substrate asperities with 
height smaller than h. 

against the nitrile rubber surface, lubricated by vari- 
ous types of oil with kinetic viscosities ranging from 
u = 5.5 X 10~^ to 46 X 10~^ m^/s. The nominal squeezing 
pressure p ^ 1 MPa. The steel cylinder was oscillating 
along its axis with the frequency / = 50 Hz and the am- 
plitude a = 1 mm. Thus, the average sliding velocity 
V ^ (jjal^j2 ^ 0.2 m/s (where uj = 27r/ is the angu- 
lar frequency). The steel surface had the rms roughness 
1.03 jam. The experiment showed that the kinetic fric- 
tion /ik ^ 0-6 was nearly independent of the lubrication 
oil viscosity. This indicate that the lubrication oils are 
nearly completely removed from the steel-rubber asperity 
contact regions and that the observed friction is due to 
the asperity induced deformation of the rubber surface 0. 
For perfectly smooth surfaces, the thickness h of the oil 
film in the contact region is easy to estimate using the 
equation 25, {puA/pty/^, where o is the oil mass 

density, u the kinematic viscosity, A the nominal contact 
area, and t the squeezing time which in the study pre- 
sented in Ref. 24] was 15 minutes. In the present case 
this gives h ^ 20 — 60 nm. The actual film thickness at 
the asperity contact regions (see Fig. IT^ may be much 
smaller. 

We now show that the rubber friction theory devel- 
oped in Ref. 2] , where the whole friction is attributed to 
the asperity-induced viscoelastic deformations of the rub- 
ber, can in fact explain the magnitude of the friction ob- 
served by Mofidi et al[2^ (which is similar to the friction 
observed for rubber seals during start-up of sliding^). 
Since the sliding surfaces are embedded in oil we will as- 
sume that there is negligible adhesive interaction between 
the surfaces so that the area of real contact is determined 
mainly by the applied pressure. In Fig. El we show the 
calculated kinetic friction coefficient /ik as a function of 
the logarithm of the sliding velocity of the rubber block. 
In the calculation we have assumed styrene butadiene 



(SB) copolymer with 60% carbon black and T = 40 °C. 
The friction is entirely due to the surface roughness of the 
substrate which is assumed to be self affine fractal with 
the fractal dimension = 2.3 and the root-mean-square 
roughness 1 /im. The long-wavelength roll-off wavevector 
qo = 27r/Ao with Aq = 100 /im, and the short wavelength 
cut-off wavevector qi = 27r/Ai where Ai = 2 nm (top 
curve), 10 nm (middle curve) and 50 nm (bottom curve). 
The solid curves are the full calculation while the dashed 
curves are without the flash temperature effect 0, E^ . 
For the Ai = 50 nm case, according to the calculation 
the area of real contact between the rubber and the sub- 
strate is only of order 0.3% of the nominal contact area so 
the frictional shear stress (at the shortest scale ~ 50 nm) 
is of order /ikP/0.003 ^ 100 MPa which is close to the 
(ideal) rupture stress of rubber at short length scale (the 
rupture stress at the macroscopic scale is usually much 
smaller due to "large" crack-like defects in most rubber 
objects). Thus, one may expect rubber wear to occur 
when sliding under the present conditions, as was also 
observed in the study reported on in Ref. 24] . 

In the context of rubber sealing, the surfaces of the 
hard countermaterials (usually steel) are prepared to 
have surface roughness with a rms amplitude in the range 
^0.1 — 1 /im, which will trap lubrication oil at the in- 
terface. After a long time of stationary contact lubricant 
fluid will only occur in the substrate valleys, but at the 
onset of sliding lubrication fluid is dragged out from the 
valleys (or cavities) to form (at high enough sliding veloc- 
ity) a thin lubrication film between the substrate asper- 
ities and the rubber surface; this will reduce the friction 
and the rubber wear. 

The theories presented above predict that the rubber 
friction decreases for high enough sliding velocity. This 
may result in stick-slip motion (as often observed, e.g., 
for rubber wiper blades), and is a major problem in many 
practical applications as it may generate strong noise and 
wear, and may lead to device malfunction. 

We note, as a curiosity, that the increase in the tem- 
perature in the asperity contact regions (flash tempera- 
ture) increases the friction at high sliding velocity, see 
Fig. El This is opposite to some other rubber appli- 
cations, e.g. for tires, where the flash temperature de- 
creases the friction 0, This difference is due to the 
fact that in the present case the temperature effect be- 
comes important only at relative high sliding velocities, 
where the asperity-induced perturbing frequencies occur 
on the glassy-side of the viscoelastic loss peak (the tan(5- 
peak); as the temperature increases the rubber viscoelas- 
tic spectra shift towards lower frequencies, i.e., the per- 
turbing frequencies will occur closer to the maximum of 
the tan(5-curve, and the friction will increase. 

We note that many earlier studies of the adhesional 
contribution to rubber friction have used "polished" 
(steel) surfaces j2lj. But such surfaces typically have sur- 
face roughness on many different length scales with a 
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FIG. 16: The nominal frictional shear stress (Ff/Ao) (a) and 
the projected (on the xy-plane) rubber-substrate contact area 
A (normalized by the nominal contact area Aq) (b) during 
sliding of a rubber block on substrates with increasing root- 
mean-square roughness /irms- (Schematic.) 

rms amplitude of order 1 /im and cannot be considered 
as smooth with respect to rubber friction. In fact, as dis- 
cussed above, the dominant contribution to the friction 
in these cases may be derived from the substrate asperity 
induced (viscoelastic) deformation of the rubber surface. 

Finally, let us discuss (qualitatively) the magnitude of 
rubber friction as a function of increasing surface rough- 
ness. At small surface roughness, because of the rubber- 
substrate adhesive interaction, there will be complete 
contact between the rubber and the substrate and the 
rubber friction will be the sum of the contributions il- 
lustrated in Fig. 2 (a) and (c), and the friction will 
increase with the surface roughness amplitude. As the 
roughness increases above some critical value, incomplete 
contact will occur at the sliding interface and when the 
surface roughness becomes large enough the friction may 
decrease. However at large enough roughness the rubber 
friction may increase again, at least if the (perpendicular) 
squeezing pressure is high enough. It is clear that even 
in the absence of a lubricant, for sealing- applications the 
best substrate surface may have some (small) roughness 
in order to minimize the friction. 

7. Summary and conclusion 

We have studied the sliding friction for viscoelastic 
solids, e.g., rubber, on hard flat substrate surfaces. We 
have shown that the fluctuating shear stress, which re- 
sult from the thermal motion of the atoms or molecules 



in a viscoelastic solid, gives rise to very strong stress- 
fluctuations, which at the nanoscale are in the MPa- 
range. This is similar to the depinning stresses which 
typically occur at solid-rubber interfaces, indicating the 
crucial importance of thermal fluctuations for rubber fric- 
tion on smooth surfaces. We have developed a detailed 
model which takes into account the influence of thermal 
fluctuations on the depinning of small contact patches 
(stress domains) at the rubber-substrate interface. The 
theory predict that the velocity dependence of the macro- 
scopic shear stress has a bell-shaped form, and that the 
low-velocity side exhibit the same temperature depen- 
dence as the bulk viscoelastic modulus, in qualitative 
agreement with experimental data. Finally, we have dis- 
cuss the influence of small-amplitude substrate roughness 
on rubber sliding friction and shown that in typical appli- 
cations to rubber sealing, the substrate asperity-induced 
viscoelastic deformations of the rubber surface may give 
the dominant contribution to the sliding friction. 
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Appendix A: Memory friction and fluctuating 
force 

Eq. (11) is a standard result in the general theory of 
Brownian motion which can be derived in various ways. 
The most general proof is based on the memory func- 
tion formalism as described, e.g., in the book by D. 
Forster 29j. A simpler (but less general) derivation of Eq. 
(11) involves the study of a particle (coordinate q{t)) cou- 
pled to an infinite set of harmonic oscillators (the heat 
bath) (coordinates x^). For the readers convenience, we 
briefly review this derivation here. The particle and the 
heat bath coordinates satisfies the equations of motion 
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and similar for q and we get from (A2) 

and from (A3), 

{f^{u)U{u')) = m^r]^kBTS{uo + c^O V A- (^5) 
From (A4) we get 

a^x^{uj) = -f{uj)q{uj) - f{uj) (A6) 
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= 



and 



(^7) 



(AS) 



Using (Al) and A(6) gives 

-mu;^q{u;) + -f{u;)q{u;) = f{uj) (A9) 
Using (A8), (A7) and (A5) it is easy to show that 



Im7(w)5(w+w') (AlO) 



Note also that (A9) is equivalent to 



mq+ dt' ^{t-t')q{t') = fit) 

J — 00 

In Sec. 3.1 we studied a system of many coupled dynam- 
ical variables 

mqi + / dt' jij{t - t')qj{t') = h{t) 



but this problem can be reduced to the problem stud- 
ied above by forming new dynamical variables, as linear 
combination of the old dynamical variables qi^ chosen so 
that becomes diagonal. 
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